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FOREWORD 


The  contract  was  initiated  under  Project  No.  1369,  "Mechanical 
Subsystems , "  Task  No.  136908,  "Landing  Gear  And  Alighting  Shock 
Absorption  Systems  And  Aircraft  Flotation."  Accomplishment  of  this 
effort  was  made  possible  wholly  through  the  use  of  Laboratory 
Director's  Fund.,  This  report  was  prepared  by  the  Aerospace  Mechanics 
Group  of  the  University  of  Dayton  Research  Institute  under  USAF 
Contract  AF  33(615)-3199.  The  work  was  conducted  under  the  direction 
of  the  Vehicle  Equipment  Division,  Air  Force  Flight  Dynamics  Labora¬ 
tory,  Wright-Patterson  Air  Force  Base,  Ohio,  Mr.  George  Sperry  (FDFM) , 
Project  Engineer. 

This  report  covers  work  conducted  from  15  May  1967  to  15  May 

1968. 

The  author  wishes  to  thank  Mr.  Sperry  for  providing  considerable 
insight  and  information  related  to  specific  Air  Force  requirements 
for  forward  area  aircraft  operations,  and  to  Mr.  Edward  Pekarek  for 
his  assistance  in  programming  the  dynamic  response  solution.  This 
report  was  submitted  by  the  author  in  May  1968. 

Publication  of  this  technical  report  does  not  constitute  Air 
Force  approval  of  the  report's  findings  or  conclusions.  It  is 
published  only  for  the  exchange  and  stimulation  of  ideas. 


A I VARS  V.  PETERSONS 

Acting  Chief,  Mechanical  Branch 

Vehicle  Equipment  Division 

Air  Force  Flight  Dynamics  Laboratory 
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ABSTRACT 


The  determination  of  aircraft  flotation  and  operation  capability  on 
semi- and  unprepared  soil  runways  is  a  critical  factor  in  developing  forward 
area  airfields  in  limited  warfare  situations.  An  investigation  was  conducted 
to  determine  the  variables  which  significantly  influence  aircraft  performance 
when  operating  on  soil  runways. 

Analysis  of  available  experimental  drag-sinkage-velocity  data  led  to 
the  defining  of  at  least  three  distinct  regions  for  which  the  sinkage  ratio- 
velobity  relationship  shows  a  distinct  response.  A  drag  ratio-sinkage  ratio 
equation  was  developed  for  use  in  one  of  these  velocity  regions.  The  effects 
of  twin  wheel  and  tandem  wheel  arrangements  were  analyzed  on  a  preliminary 
basis.  The  results  of  a  sinkage  study  using  presently  available  sinkage  pre¬ 
diction  equations  indicated  that  present  sinkage  analysis  accuracy  is  in  the 
range  of  ±50%  to  ±100%. 

In  order  to  develop  a  suitable  flotation  criteria,  an  investigation  was 
conducted  of  the  dynamic  landing  gear  contacting  element-soil  interaction 
response  utilizing  elastic  theory.  These  results  led  to  the  development  of  a 
flotation  parameter  (related  to  sinkage)  and  a  flotation  index  (related  to  drag) 
in  nomographic  form  which  permits  comparative  flotation  analysis  of  landing 
gear  systems. 

This  document  is  subject  to  special  export  controls  and  each  transmittal 
to  foreign  governments  or  foreign  nationals  may  be  made  only  with  prior 
approval  of  the  Air  Force  Flight  Dynamics  Laboratory  (FDFM),  Wright- 
Patterson  Air  Force  Base,  Ohio  45433.  The  reason  for  this  distribution 
limitation  is  that  the  report  would  enhance  Sino-Soviet  bloc  military  and 
industrial  research  programs  having  a  potential  strategic  value. 
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SECTION  I 
INTRODUCTION 

The  conduction  of  limited  warfare  in  many  cases  requires  the 
operation  of  aircraft  on  forward  area  semi-  and  unprepared  soil  surfaces. 
The  effective  and  efficient  utilization  of  aircraft  demands  that  landings 
on,  maneuvers  on,  and  takeoff  from  minimum  length  and  width  soil 
surfaces  be  accomplished  with  maximum  aircraft  control  and  safety  at 
maximum  loads  with  the  smallest  possible  deterioration  and  damage  to 
the  aircraft  and  airfield  surface. 

Since  of  necessity  these  forward  airfields  are  limited  in  ground 
strength  and  runway  length,  aircraft  landing  gear  sinkage,  ground  drag, 
and  surface  roughness  are  critical  factors  to  the  utilization  of  soil  sur¬ 
faced  airfields.  An  indication  of  the  magnitudes  of  sinkage  is  shown  in 
Figure  1.  Currently  available  flotation  criteria  and  analysis  is  based  on 
empirical  results  of  full,  scale  experimental  tests  conducted  with  only 
limited  ranges  of  such  variables  as  load,  tire  parameters,  landing  gear 
configuration,  and  soil  conditions. 

The  basic  objective  of  this  program  is  to  analytically  define  landing 
gear  -  soil  surface  interaction  in  order  to  establish  adequate  aircraft 
flotation  criteria.  This  objective  is  to  be  accomplished  by  a  three  phase 
program.  Phase  I,  described  in  this  report,  included  a  determination  of 
the  critical  parameters  affecting  flotation,  an  investigation  of  available 
experimental  data  and  analyses  for  defining  landing  gear  sinkage,  surface 
drag,  and  the  interrelationship  between  sinkage  and  drag  leading  to  the 
ievelopment  of  a  suitable  mathematical  model  for  establishing  flotation 
criteria. 

Phar  e  II  is  concerned  with  the  development  of  analysis  equations 
for  use  with  computer  simulation  techniques  to  analytically  define  the  range 
of  variables  associated  with  flotation.  Initial  development  will  be  made  of 
a  sinkage  and  drag  index  for  defining  aircraft  flotation  in  nomographic 
form.  In  Phase  III  the  effects  of  braking,  high  velocity,  and  multiple 
wheel  effects  will  be  analyzed  leading  to  the  development  of  a  final  flotation 
criteria  for  defining  aircraft  performance  on  semi-  and  unprepared  soil 
airstrips. 


Figure  1  Aircraft  Tire  Sinkage 


SECTION  II 


IDENTIFICATION  OF  VARIABLES 

1.  Literature  Review 

In  studies  conducted  to  date  on  mobility  and  flotation,  numerous 
variables  have  been  postulated  as  influencing  aircraft  response  and  landing 
gear  -  soil  interaction  for  aircraft  operating  on  semi-  and  unprepared  soil 
runways.  Existing  literature  discusses  these  variables  in  terms  of  traction 
wheels  and  transport  wheels.  Since  the  behavior  of  an  aircraft  tire  is 
essentially  that  of  a  transport  wheel,  only  those  variables  associated  with 
the  transport  wheel  will  be  treated.  Freitag  (1)  in  a  comprehensive  review 
of  available  mobility  data,  included  tire  width,  tire  diameter,  tire  inflation 
pressure,  and  the  cone  index  (soil  parameter)  as  constituting  the  primary 
variables  related  to  drawbar  pull  (drag)  in  single  wheels.  Nuttal  (2)  included 
the  additional  variables  of  load,  tire  deflection,  vehicle  speed,  soil  viscosity, 
soil  cohesion,  and  soil  friction  in  analyzing  drag  and  sinkage  of  wheels. 

Other  researchers  (3,4,5)  considered  essentially  the  same  array  of  variables 
in  analyzing  mobility. 

2.  Analysis  of  variables 

The  numerous  variables  associated  with  the  landing  gear  -  soil  interaction 
all  to  some  extent  influence  the  mechanisms  of  response  in  dynamic  load 
situations.  Looking  first  at  the  characteristics  of  soil  (ground)  and  the 
landing  gear  tires,  Table  I  gives  a  list  of  the  significant  variables  commonly 
specified,  without  at  this  stage  involving  any  interaction  phenomena. 

TABLE  I 

Significant  Variables:  Soil,  Landing  Gear  Tires 

Landing  Gear  Tires 


Contact  Area  and  Geometry* 

Inflation  Pressure 
Tire  Ply  Rating 
Tire  Stiffness 
Tire  Deflection 
Contact  Stress  Distribution 
Diameter  and  Width  of  Tire 
Number  of  Tires 
Tire  Spacing  and  Configuration 
Tread  Pattern 
Coefficient  of  Friction:  Tire-Ground 

*  Other  landing  gear  variables  associated  with  landing  gear  -  soil  interaction 
not  considered  in  this  study  include:  landing  gear  stiffness,  braking 
system,  antiskid  and  steering  system  effects. 


Soil 

Elastic  Failure 

G=  Shear  Modulus  cs=  Cohesion 
p=  Poisson's  Ratio  <p=  Friction 
p-  Mass  Density 
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Upon  introducing  dynamic  loading  conditions  and  relating  the  problem 
to  interaction  phenomena,  the  number  of  variables  increases  significantly. 
Figure  2  provides  a  systematic  overall  look  at  the  interaction  problem. 
Analysis  of  e>dsting  literature  indicates  that  the  suitability  of  soil  airstrips 
for  aircraft  operations  is  directly  related  to  vertical  sinkage  of  the  air¬ 
craft  wheel,  ground  drag,  and  roughness  of  the  runway, surface.  On  this 
basis,  the  interaction  problem  can  be  defined  in  terms  of: 

(1)  Loads 

(2)  Soil  Characteristics 

(3)  Vehicle  Characteristics 

(4)  Surface  Characteristics 

The  most  important  property  of  soil  is  its  deformability  characteristics 
in  loading  and  unloading.  Equations  of  state  relating  load  to  sinkage  have 
been  developed  by  a  number  of  researchers  in  mobility  (1,3,4). 

Directly  related  to  the  loading  and  unloading  characteristics  of  soil 
is  the  residual  permanent  deformation  for  each  wheel  passage.  This 
permanent  deformation,  when  combined  with  the  initial  runway  roughness, 
determines  the  ground  roughness  for  subsequent  aircraft  operations. 

In  addition  to  the  deformability  of  soil,  the  tire  characteristics  also 
influence  the  resultant  sinkage  of  an  aircraft  wheel  under  fixed  loading 
conditions.  The  numerous  variables  associated  with  tire  characteristics, 
as  indicated  in  Figure  2,  can  be  thought  of  in  terms  of  a  flotation  parameter, 
F,  related  to  sinkage.  Thi3  flotation  parameter  could  then  be  used  to 
compare  landing  gears  (all  other  conditions  being  constant  except  the  tire 
characteristics)  to  provide  a  rational  means  for  selecting  the  higher 
flotation  gear  system. 

The  third  significant  quantity  in  defining  suitability  of  a  landing  gear 
for  operation  at  a  site  is  ground  drag  on  the  aircraft  wheel  due  to  the 
wheel  -  ground  interaction.  While  drag  is  helpful  in  a  landing  operation, 
large  drag  magnitudes  in  takeoff  operations  leads  to  excessive  runway 
length  requirements.  The  possibility  also  exists  of  immobilization  of  the 
aircraft.  As  indicated  in  Figure  2,  ground  drag  is  influenced  by  load, 
soil,  and  vehicle  characteristics.  Subsequent  analysis  has  indicated  that 
over  a  significant  portion  of  the  aircraft  horizontal  velocity  range  the 
dominant  parameter  in  determining  drag  is  aircraft  wheel  sinkage. 
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Figure  Z  Variables  in  Landing  Gear-Unprepared  Runway  Problem 
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The  large  number  of  variables  listed  in  Figure  2  can  be  interpreted 
in  terms  of: 

(1)  loads 

(2)  deformability  of  soil 

(3)  flotation  parameter  and  flotation  index 

(4)  initial  roughness 

if  sinkage,  drag,  and  i*oughness  are  utilized  as  the  significant  quantities 
in  defining  the  suitability  of  an  unprepared  site  for  aircraft  operations. 


SECTION  III 

DRAG-SINKAGE  ANALYSIS 

1.  Drag  Analysis 

The  drag  forces  exerted  by  the  deformed  ground  on  the  aircraft  wheel 
are  an  important  consideration  in  relation  to  aircraft  thrust  for  determin¬ 
ing  the  takeoff  capability  and  runway  length  for  aircraft  operations  on  soil 
surfaces.  Figure  3  shows  the  pneumatic  wheel  -  ground  interaction 
phenomena  resulting  in  drag  forces  being  exerted  on  the  wheel.  In  this 
report,  sinkage,  Z,  is  defined  as  the  instantaneous  displacement  of  the 
ground  beneath  the  tire  from  its  equilibrium  position  (see  Figure  3). 
Depending  upon  the  inflation  pressure,  a  pneumatic  tire  may  act  as  either 
a  rigid  wheel  or  a  nonrigid  (flexible)  wheel  on  soft  ground.  The  numerous 
mobility  and  flotation  studies  (1,2,  3,  6)  conducted  in  the  past  decade  have 
provided  considerable  empirical  and  semi-theoretical  information  on  the 
drag  characteristics  of  rigid  and  flexible  wheels  on  rigid  surfaces.  Con¬ 
siderably  less  information  is  available  for  wheel  drag  on  soft  surfaces. 

2.  Drag  Variables 

The  pertinent  variables  affecting  drag  as  shown  in  Figure  2  are 
summarized  in  Table  II. 

TABLE  II 
Drag  Variables 

Friction 

Quantity  Drag  Load  Slip  Cohesion  Angle 
Symbol  R  P  S  c  co 

Units  F  F  FL*2 

Shear  Mass  Soil 

Modulus  Density  Viscosity 

G  P  n 

-2  -4  2  -2 

FL  FL  T  FL  T 

Contact  Tire  Tire  Coefficient 

Velocity  Acceleration  Area  Width  Diameter  of  Friction 

V  g  A  b  D  \ 

lt”1  lt"2  l2  L  L 


Impingement 
Drag  Coefficient 
C 


Poisson's 

Ratio 


The  effect  of  horizontal  aircraft  velocity  on  drag  is  indicated  in 
Figure  4  which  shows  the  experimentally  measured  trend  of  the  drag 
ratio  (Drag/Vertical  Wheel  Load)  with  horizontal  velocity  (5).  Reference 
to  Figure  4  indicates  at  least  three  distinct  zones  of  drag  ratio 
response.  At  low  speeds  (less  than  5  knots),  the  drag  ratio  changes 
considerably  with  velocity  due  to  the  viscous  effects  of  the  soil  as  related 
to  tire  sinkage.  From  5  knots  to  approximately  40  to  50  knots,  the  drag 
ratio  is  relatively  constant,  indicating  that  the  rate  of  loading  effects  in 
the  soil  are  no  longer  significant.  In  the  third  region  (greater  than 
approximately  40  to  50  knots),  the  drag  ratio  increases  rapidly  with 
increasing  velocity  and  then  diminishes  apparently  due  to  the  combined 
effects  of  soil  inertia  drag  and  hydroplaning  effects.  It  should  be  recognized 
that  even  though  the  drag  ratio  increases  rapidly  in  the  third  region,  the 
resultant  drag  force  (which  is  the  critical  factor  in  determining  runway 
length)  very  likely  remains  less  than  the  drag  in  region  two,  due  to  the 
large  lift  effects  at  the  higher  horizontal  velocities  and  consequent  reduction 
in  vertical  wheel  loads. 

Assuming  that  immobility  under  static  loads  is  not  a  problem,  the 
most  critical  drag  effects  occur  in  region  two  and  three.  Other  than 
limited  experimental  data  indicating  the  variation  of  the  drag  ratio  with 
horizontal  velocity,  little  information  is  available  for  understanding  the 
interrelationship  between  drag,  sinkage,  high  velocity,  and  soil  inertia. 

On  the  assumption  that  a  tire  traversing  soil  at  high  velocities  responds 
the  same  as  a  tire  traversing  slush  at  high  velocities,  Boeing  (5)  utilized 
the  drag  equations  developed  for  slush  conditions  to  interpret  high  velocity 
soil  drag  in  the  form: 

D'=KC  ZV8  (1) 

where  D'=  additional  drag  due  to  soil  inertia  effects 
K  =  constant  related  to  soil  density  (p)  and 
tire  width  (b)  =  1/2 pb 

C^  =  impingement  drag  coefficient  related  to  soil 

Z  =  instantaneous  sinkage 
V  =  horizontal  velocity 

The  drag  magnitude  defined  by  Equation  1  must  be  added  to  the  drag 
developed  in  region  two  to  determine  the  total  drag  acting  on  the  tire  in 
region  three.  Very  little  experimental  data  or  theoretical  justification  is 
available  for  proving  or  disproving  the  validity  of  Equation  1.  Also,  the 
soil  impingement  drag  coefficient  is  not  known  for  varying  soil  types. 

More  recently,  Lockheed  (6)  has  begun  high  speed  drag  ratio  tests  which 
should  provide  additional  information  for  establishing  a  relationship  between 
drag,  sinkage,  and  soil  inertia  at  high  velocities. 


9 


VERTICAL  LOAD 


*  ^  JTgVafeg^<Sfe^«sg^  -y.yjw*a 


I: 


£ 

I 


Figure  4  Drag  Ratio  vs  Horizontal  Velocity  on  Soil 
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3.  Drag-Sinkage,  Single  Wheel 

In  region  two,  experimental  evidence  exists  which  indicates  that  the 
effects  of  velocity,  slip,  viscosity,  soil  density,  and  soil  inertia  have 
only  a  small  influence  on  drag  characteristics.  The  roiling  resistance 
at  zero  sinkage  is  also  small  in  comparison  to  the  total  drag  at  sinkages 
normally  encountered  in  soil. 


A  dimensional  drag  analysis  of  the  remaining  variables  in  Table  II 
could  take  the  form: 


R/P=  / 


Additionally,  if  it  is  recognized  that  the  primary  factors  related  to 
sinkage,  Z,  are  the  load  (P),  geometry  factors  (b,A),  and.  the  soil  propertie 
(p,  G,  c,  cp),  then  a  simplified  drag  ratio  relationship  of  the  form  given  in 
Equation  3  is  obtained. 


R/P=/(Z/D)  (3) 

The  inclusion  of  the  wheel  diameter,  D,  in  the  drag  ratio  equation  is 
based  on  simple  work-energy  relationships  which  provide  an  approximation 
to  drag  analysis  for  region  two  velocity  rolling  pneumatic  tires.  Reference 
to  Figure  3  indicates  that  the  work  required  to  increase  the  potential  energy 
of  the  wheel  as  it  moves  from  position  1  to  position  2  is  R*  J i  where  l  is 
related  to  the  diameter,  D,  and  consequently  can  be  written  as:  R  =  KiD. 

The  soil,  however,  undergoes  a  deformation,  the  sinkage  Z  (see  position  3) 
which  results  in  the  work  quantity  (or  energy  loss):  KgPZ.  Neglecting  any 
other  forms  of  energy,  the  work-energy  equation  becomes 


RK1D=I^PZ 

(4) 

or 

R/P=  /(Z/D) 

(5) 

Such  a  simplified  drag  ratio  relationship  can  also  be  developed  by 
reference  to  existing  mobility  literature.  Bekker  (3)  has  examined  the 
relationship  among  load,  tire  diameter,  and  drag  for  pneumatic  tires  on 
soil  using  a  load-sinkage  relationship  of  the  form 

i  \ 

p=  kc  +  k  J  Zn 

U  j 

(6) 

where  p  =  contact  stress 

k  ,  k  =  soil  deformation  moduli 
c  cp 

n  =  exponent  of  sinkage 

b  =  smallest  dimension  of  contact  area 
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For  the  case  of  n  =  1,  and  neglecting  the  rolling  resistance  at  zero  sinkage, 
the  drag  equation  as  given  by  Bekker  (3)  becomes 


where 


b  (pn-y 

2  <kc  +  by 


p.  =  inflation  pressure 
p^  =  carcass  pressure 


The  vertical  equilibrium  force  on  the  tire  is  approximately  given 
P=b^1{pj.+pc)  (8; 


where 


4=  contact  length  between  tire  and  ground 
P  =  total  load  on  tire 

Combining  equations  (6)  and  (7)  and  noting  that 


leads  to 


Mp.  +  pc) 

Z  =  (k  +  bk  ) 

C  CD 


R/P  = 


Recognizing  that  is  related  to  the  diameter  of  the  tire,  D,  (i.  e. ,  for 
constant  tire  deflection,  varies  with  D)  then  Equation  10  can  be  written  as 

R/P  =  /(Z/D)  (11) 

While  it  might  be  reasoned  that  Xj  is  related  to  both  the  tire  diameter  and 
the  tire  deflection,  d,  the  effects  of  the  tire  deflection  are  inherent  in  the 
sinkage,  Z,  which  includes  the  inflation  pressure,  p^.  The  form  of  the  drag 

ratio  relationship  as  given  in  Equation  11  can  also  be  established  from  a 
mobility  analysis  by  Nuttal  (2). 

On  the  basis  of  the  above  analysis  and  a  review  of  the  literature,  an 
investigation  was  made  of  the  relationship  between  drag  and  sinkage  based 
on  existing  experimental  evidence.  While  numerous  drag  studies  have  been 
conducted,  the  measurement  of  sinkage  in  many  cases  was  not  taken. 

Figure  5  presents  the  drag  ratio  (R/P)  versus  sinkage  ratio  (Z/D)  based  on 
data  from  four  sources  (1,2, 3,  7).  The  plotted  data  includes  the  full  range 
of  soil  type  (sand  to  clay)  and  tire  diameters  from  approximately  30”  to  70" 
which  encompasses  the  range  of  tire  sizes  for  heavier  military  aircraft.  In 
reference  to  Figure  5,  it  is  clearly  evident  that  there  is  an  increasing  trend 
of  the  drag  ratio  with  increasing  sinkage  ratio  for  Z/D  varying  from  0  to  C.  10. 
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Little  interest  is  given  to  the  relationship  beyond  a  Z/D  of  0.  10,  since  for  a 
30"  diameter  tire,  the  sinkage  would  be  in  excess  of  3"  which  is  presently 
the  instantaneous  sinkage  normally  specified  as  a  failure  sinkage  for  military 
aircraft  (9). 

It  is  also  evident  from  analysis  of  Figure  5  that  there  is  a  considerable 
amount  of  scatter  in  the  data.  This  scatter  may  arise  from  several  causes. 
The  full  range  of  soil  type  is  included  in  the  analysis.  As  additional  drag- 
sinkage  data  becomes  available,  the  effect  of  soil  type  may  be  more  clearly 
recognized.  In  some  instances  the  data  sources  do  not  specify  whether  the 
definition  of  sinkage  does  or  does  not  include  the  bow  wave,  H  (see  Figure  3). 
Sinkage  as  utilized  in  this  report  does  not  include  the  bow  wave  (H). 
Additionally,  soil  viscosity  effects  may  be  present  to  some  extent  for  those 
test  data  obtained  at  the  lower  velocities.  Also  it  is  possible  that  the  sinkage 
quantity,  Z,  may  not  adequately  reflect  the  variables  b  (tire  width),  and  d 
(tire  deflection). 

In  general,  however.  Figure  5  suggests  that  the  most  significant 
variables  as  related  to  drag  are  the  load  (P),  sinkage  (Z),  and  tire  diameter 
(D).  Consequently,  a  least  squares  linear  and  a  least  squares  parabola  were 
fitted  to  the  data  with  the  sinkage  ratio  as  the  independent  variable  and  the 
drag  ratio  as  the  dependent  variable.  Analysis  of  the  results  indicated  that 
the  quadratic  term  was  not  a  significant  contributant  to  the  drag  relationship 
in  the  range  of  Z/D  considered  and  therefore  based  on  the  experimental  data, 
the  drag  ratio  -  sinkage  ratio  relationship  is  given  by 

R/P  =  0.015  +  3.09  (Z/D)  (12) 

The  positive  R/P  intercept  at  zero  sinkage  would  be  expected  since  some 
rolling  resistance  (friction)  is  present  at  zero  sinkage. 

Of  additional  interest  is  the  determination  of  the  upper  bound  drag 
magnitudes  likely  to  be  encountered  by  an  aircraft  on  an  unprepared  landing 
site  when  the  characteristics  of  the  aircraft  and  approximate  sinkage  values 
(based  on  soil  type  and  conditions)  are  known.  For  this  purpose  80%  and  90% 
probability  lines  were  established  and  are  shown  in  Figure  5.  This  upper 
bound  drag,  when  related  to  aircraft  thrust  capability,  should  be  helpful  in 
determining  runway  lengths  and  possible  immobility  conditions  for  aircraft 
operating  on  semi-  and  unprepared  land  sites. 

Subsequent  to  the  least  squares  analysis  used  to  establish  Equation  12, 
additional  drag  data  became  available  from  tests  reported  by  Boeing  (5).  This 
additional  data  is  shown  in  Figure  5  and  tends  to  verify  the  established 
relationship. 

4.  Drag-Sinkage,  Twin,  Tandem  Wheels 

Very  little  information  (theoretical  or  experimental)  is  available 
concerning  multiple  wheel  effects  as  related  to  drag  forces  on  aircraft 
wheels.  Two  types  of  multiple  wheel  configurations  are  of  interest:  dual 


and  tandem.  Dual  wheels  at  very  large  spacings  tend  to  act  as  single 
isolated  wheels.  The  drag  relationship  for  closely  spaced  wheels  is  not 
Known.  Similarly,  the  drag  relationship  for  tandem  wheels  where  one 
wheel  follows  in  the  path  of  the  first  wheel  has  not  been  adequately  determined. 

In  analyzing  multiple  wheel  effects  on  drag  based  on  existing  experi¬ 
mental  data  (1965),  Freitag  (1)  concluded  that  in  frictional  soils  existing 
data  indicated  that  each  wheel  in  a  dual  pair  performs  at  the  same  level  as 
an  isolated  single  wheel  for  all  dual  wheel  spacings.  Analysis  of  dual  wheel 
effects  in  frictionless  soils  permitted  no  conclusions  due  to  insufficient  and 
inconsistent  data.  Even  the  conclusions  drawn  from  studies  in  frictional 
soils  are  questionable,  due  to  the  limited  test  data  and  the  fact  that  the  test 
data  covered  only  two  tires  with  a  range  in  wheel  diameter  from  24"  to  40", 
while  dual  tire  spacing  ranged  from  8"  to  16". 

A  limited  amount  of  additional  test  data  involving  drag  measure¬ 
ments  for  multiple  wheels  is  available  from  a  recent  WES  (7,  9)  flotation 
investigation  on  clay  soil.  Twin  wheel  and  multiple  wheel  (twin-twin  or 
twin-tandem)  configuration  tests  were  conducted.  An  indication  then  of 
multiple  wheel  effects  can  be  gained  by  analyzing  the  drag  sinkage  data  from 
only  the  WES  study.  For  comparative  purposes  a  linear  least  squares  fit  of 
the  WES  single  wheel  data  yields  the  equation  (see  Figure  6) 

R/P  =0.030  +  1.65  (Z/D)  (13) 

Due  to  the  limited  amount  of  twin  wheel  test  data,  all  the  twin  wheel  drag- 
sinkage  data  has  been  plotted  in  Figure  7  on  the  basis  of  taking  the  total  drag 
and  dividing  by  the  total  of  the  wheel  loads  (converting  R/P  to  a  single  wheel 
basis).  The  twin  wheel  spacing  ranged  from  25"  to  60".  A  linear  least 
squares  analysis  of  this  test  data  yields 

R/P  =  0.  030  +  1.  45  (Z/D)  (14) 

There  is  no  apparent  trend  in  the  larger  spaced  wheels  data  of  lying  above 
the  linear  least  squares  line  as  would  be  expected  if  the  wheels  tended  to  act 
as  single  isolated  wheels  at  larger  spacings.  Comparison  of  Equations  13 
and  14  (single  wheel  to  twin  wheel  respectively)  indicates  that  while  there  is 
some  effect  of  wheel  spacing  o:,i  the  drag  ratio-sinkage  ratio  relationship 
(less  than  15%  in  the  range  considered),  it  is  not  a  pronounced  effect.  Since 
most  conventional  aircraft  tires  for  cargo  aircraft  cannot  be  spaced  closer 
than  about  20  inches,  it  appears  that  each  wheel  in  a  twin  wheel  configuration 
tends  to  act  as  a  single  isolated  wheel  in  frictionless  soils  in  determining 
drag  characteristics. 

Multiple  wheel  (three  or  more  wheels  in  both  tandem  and  dual 
configurations)  drag-sinkage  data  taken  from  the  WES  tests  is  shown  in 
Figure  8.  A  linear  least  squares  analysis  of  this  data  gives 

R/P  =  0.  020  +  1,  26  (Z/D)  (15) 


15 


There  is  a  considerable  difference  between  Equations  15  (multiple  wheel) 
and  Equation  13  (single  wheel).  This  difference  is  due  to  tandem  effects; 
that  is,  wheels  following  in  the  same  path  as  forward  wheels.  Since  the 
lead  tire  causes  ground  deflection  (with  subsequent  changes  in  soil  deform- 
ability  properties)  subsequent  to  the  follow-on  tire,  the  tandem  tire  will 
undergo  less  sinkage  and  consequently  less  drag. 

Based  on  the  limited  data  presented  in  Figures  6,  7,  and  8,  it  is 
difficult  to  conclude  an  exact  drag-sinkage  relationship  in  frictionless  soils 
for  twin  and  tandem  wheel  configurations.  It  would  appear,  however,  that 
for  normal  tire  spacings,  dual  wheels  tend  to  act  as  single  isolated  wheels 
in  determining  drag,  whereas  a  significant  reduction  in  drag  (over  that 
determined  by  considering  the  wheels  as  single  isolated  wheels)  can  be  ex¬ 
pected  for  tandem  wheel  configurations  in  frictionless  soils. 

All  of  the  data  used  in  developing  Figures  5,  6,  7,  and  8  and  the 
resulting  analysis  equations  are  summarized  in  Appendix  I. 

5.  Summary 

Based  on  the  above  analysis  and  evaluation  of  existing  experimental 
drag  data,  the  critical  parameter  in  determining  drag  on  an  aircraft  wheel 
in  the  Region  two  velocity  range  is  sinkage  of  the  tire  into  the  soil.  If  an 
estimate  can  be  made  of  the  anticipated  sinkage,  then  the  total  drag  can  be 
determined  within  suitable  limits  for  developing  required  runway  lengths  as 
related  to  specific  aircrafts.  As  additional  drag  data  becomes  available, 
the  drag  ratio-sinkage  ratio  as  defined  in  Figure  5  and  Equation  12  can  be 
suitably  modified  to  reflect  this  additional  drag  information. 

Analysis  of  aircraft  drag  ratio  response  at  high  velocities  (Region 
three)  has  not  been  adequately  defined  either  theoretically  or  experimentally 
to  permit  reliable  estimates  of  drag  magnitude.  If  the  high  velocity  drag 
relationship  as  developed  by  Boeing  (5)  is  utilized,  then  the  drag  ratio  defining 
aircraft  response  in  Region  three  can  be  treated  as  the  sum  of  the  drag  de¬ 
fined  by  Region  two  plus  the  drag  determined  from  Equation  1. 


SECTION  IV 


SINKAGE  PREDICTION  EQUATIONS 


I.  General 

Since  vertical  sinkage  of  the  aircraft  tire  is  the  dominant  factor  in 
determining  drag  on  the  wheel  in  the  region  two  velocity  range,  and  based 
on  the  work  of  Boeing  (5)  forms  the  basis  for  which  inertia  drag  effects 
must  be  added  for  defining  drag  in  the  region  three  velocity  range,  it  is 
important  to  analyze  the  currently  used  sinkage  prediction  equations  for  the 
reliability  and  accuracy  of  sinkage  prediction  as  related  to  the  loading, 
landing  gear,  and  soil  variables.  Additionally  the  sinkage  affects 
aircraft  capability  to  operate  from  a  site  due  to  resulting  permanent 
deformation  and  roughness  effects,  and  thus  forms  the  basis  of  flotation 
analysis. 

Numerous  experimental  and  analytical  efforts  have  been  undertaken 
to  develop  load-sinkage  equations  which  adequately  interpret  soil  properties 
as  related  to  soil  deformation  characteristics  under  loading.  Defining  soil 
load-sinkage  behavior  is  an  extremely  complex  problem  due  to  the  inherent 
nonlinear  and  inelastic  behavior  of  soil.  Under  low  stress  loadings,  the 
soil  response  is  approximately  elastic.  At  higher  stress  levels,  which  is 
common  to  problems  in  flotation  and  mobility,  the  contact  stresses  are  far 
in  excess  of  any  elastic  limit  for  soil.  For  aircraft  wheels  moving  over 
soft  soils,  the  complex  soil  response  involves  elastic,  viscous  (time 
dependent),  and  plastic  deformations. 


2.  Currently  Used  Sinkage  Equations 


Initial  efforts  at  developing  load-sinkage  equations  involved  power 
functions  of  the  sinkage  term,  Z.  The  most  widely  used  equation  in  this 
form  was  developed  by  Bekker  (10)  as 


where 


(16) 


p  =  contact  stress,  psi 

k  ,  k  =  soil  deformation  moduli  related  to  friction  and 
C  ^  cohesion  of  soil 
n  =  exponent  of  sinkage 

b  =  smallest  dimension  of  contact  area,  in. 

Z  =  vertical  sinkage,  in. 


Equation  16  utilizes  two  soil  parameters  to  account  for  the  frictional  and 
cohesive  contributions  to  soil  strength.  The  Bekker  developed  load-sinkage 
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equation  is  used  primarily  by  the  .Land  Locomotion  Laboratory  in  their 
studies  of  both  tire  and  track  mobility  over  soil  terrain.  Widespread  use 
of  Equation  16  has  not  developed  due  to  the  dimensional  inconsistency  in  Z 
being  a  power  function,  and  the  lack  of  a  consistent  trend  in  and  k  with 

soil  type  (negative  values  of  the  soil  deformation  moduli  have  been  reported 
in  some  instances). 

WES(5)  has  also  developed  a  load-sinkage  relationship  specifically 
for  pneumatic  tires  which  is  defined  in  terms  of  a  clay  mobility  number 
(CMN)  and  a  sand  mobility  number  (SMN).  These  mobility  numbers  as 
defined  in  Equations  17  and  18  are  developed  from  cone  index  readings 
which  are  a  measure  of  the  resistance  of  the  soil  to  penetration  by  a  cone 
tipped  rod. 

/  6  \1/S 
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Clay  Mobility  Number  (CMN)  = 


Sand  Mobility  Number  (SMN)  = 


where 

Cl  =  average  cone  index  over  first  six  inches,  psi 

G  =  slope  of  cone  index  versus  depth  averaged  over  a  depth 

equal  to  the  tire  width,  psi  per  in. 

b  =  tire  width 

=  tire  deflection,  in. 
t 

h^  =  tire  section  height,  in. 

Approximations  to  the  mean  sinkage  ratio-mobility  number  (5)  yields 

Z  —  2. 

-  =  0.  003  +  (CMN)  *  for  3  <:  CMN  £  10 


=  0.  003  +  (SMN)"1’  5  for  5  £  SMN  £  40 


Equations  19  and  20  which  define  sinkage  in  the  extreme  soil  types  (sand 
and  clay)  were  developed  by  a  dimensional  analysis  of  the  experimental  data 
taken  from  pneumatic  tire  tests  on  only  two  soils,  Yuma  sand  and  Buckshot 
clay.  Consequently  the  extension  of  Equation  19  and  20  to  other  soils  and 
conditions  is  questionable  at  the  present  time. 

An  alternate  approach  to  defining  a  load-sinkage  relationship  has 
been  developed  by  Assur  (11)  who  considered  three  possible  modes  of  soil 
response  as  given  by: 
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p  =  K  Z  1  - 
s 


(fluidization) 


(21) 


where 


Ik 

3 

(K  ) 
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1  P«  J 

Pm  j 

*  *  *  * 
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P  = 


s  •  Z 


1  -  Z3/Zs 
m 


(compaction) 


K 


P  = 


s-  Z 


1  +  Z3/Z3 
m 


(collapse) 


K  =  coefficient  of  3ubgrade  reaction,  psi  per  in. 
s 

Pm=  maximum  bearing  strength,  psi 
Z^  =  sinkage  at  maximum  bearing  strength,  in. 


(22) 


(23) 


While  these  load-sinkage  equations  can  be  used  to  approximate  almost  any 

observed  load-sinkage  response  of  soil  (through  curve  fitting  techniques),  it 

is  not  always  evident  in  attempting  to  predict  sinkages,  the  value  of  the 

maximum  bearing  strength  (p  )  and  the  sinkage-  at  the  maximum  bearing 

m 

strength  (Z  ).  Also  these  equations  are  not  related  to  any  fundamental  soil 
m 

property  which  hinders  their  utility  of  use. 

More  recently  Boeing  (5)  developed  a  load-sinkage  equation  analogous 
in  form  to  the  equation  of  Bekker.  The  nonlinear  relationship  between  load 
and  sinkage  is  accounted  for  in  Boeing's  equation  by  use  of  a  nonlinear 
instantaneous  dynamic  soil  spring  rate  (k^)  as  shown  in  Equation  24. 


where 


and 


By  letting  n  =  1  in  Equation  16,  and  defining  the  dynamic  soil  spring  rate  as 


Z  = 


Ak. 


(24) 


k^  =  kg  D  ?  =  instantaneous  dynamic  soil  spring  rate, 
^  psi  per  in. 


k  =  static  soil  spring  rate,  psi  per  in. 
s 

D  =  dynamic  factor  related  to  duration  of  loading  (taken  as 
^  unity  for  frictional  soils) 

Y  =  attenuation  factor  which  varies  with  sinkage 
A  =  contact  area,  in.  3 
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it  is  evident  that  Equations  16  and  24  are  similar  in  form. 

Because  of  the  broad  range  and  number  of  variables  selected  in 
developing  these  sinkage  relations  and  due  to  the  fact  that  the  expressions 
for  sinkage  vary  from  analytical  with  little  experimental  justification  to 
strictly  empirical  formulations,  it  is  difficult  to  define  which  of  these 
sinkage  prediction  equations  is  most  suitable  and  accurate  for  predicting 
sinkages  of  aircraft  tires  into  soil.  If  a  suitable  correlation  existed  between 
the  numerous  soil  parameters,  a  more  comprehensive  and  meaningful 
study  of  flotation  could  be  made  utilizing  experimental  data  from  many 
sources.  No  valid  relationship  has  been  established  at  the  present  time 
between  the  soil  parameters  defined  in  Equations  16,  17,  18,  21,  22,  23 
and  24. 


3.  Limited  Correlation  of  Sinkage  Prediction  Equations 

An  extension  of  the  recent  work  by  janosi  (12)  when  combined  with 

experimentally  established  relationships  between  the  California  Bearing 

Ratio  (CBR)  and  the  cone  index  (Cl)  permits  a  limited  correlation  between 

Bekker's  moduli  of  soil  deformation  and  WES's  cone  index  for  sands  and 

clays.  Based  on  the  definition  of  the  static  soi]  spring  rate,  this  correlation 

can  also  be  extended  to  include  Boeing's  soil  parameter,  k  . 

s 

The  correlation  between  soil  parameters  is  based  on  an  analysis  of 
the  cone  penetration  test.  The  cone  penetration  test  as  developed  by  WES 
(1)  consists  of  measuring  the  pressure  required  to  force  a  standard  cone 
tipped  rod  into  soil.  The  cone  index  (Cl)  then  is  defined  as  the  measured 
pressure.  The  equivalence  between  Bekker's  soil  parameters  and  WES's 
cone  index  as  developed  by  Janosi  (12)  is  given  by 


Cl  =  1.  625  s 


c 

n+1 


[(Z+1.5)n+1]-Zn+1 


40.  5175  k 


cp 


(Z+1.5) 


,n+2 


(n+l)(n+2) 


n+2 

6 

n+2 


(Z+I.  5) 
n+1 


,n+fl 


(25) 


a)  Sinkage  in  Sand 

For  sands,  k^  is  zero  and  experimental  evidence  indicates  that 
n=l  is  approximately  valid.  For  these  conditions  Equation  25  reduces  to 

Cl  =  k  (0.  947Z  +  0.  474)  (26) 
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As  indicated  in  Equations  18  and  20,  the  sand  mobility  number  (SMN)  is 
related  to  sinkage.  The  quantity  G  which  is  defined  as  the  slope  of  the  cone 
index  versus  depth  curve  averaged  over  a  depth  equal  to  the  tire  width  is 
used  to  define  the  properties  of  sand.  By  differentiating  Equation  26,  the 
quantity  G  becomes 

G  =  =  0.  947  k  (27) 

dZ  cp 

By  definition  the  cone  index  is  the  pressure  required  to  penetrate  the  cone 
tipped  rod  into  soil.  Boeing's  static  soil  spring  rate  can  then  be  developed 
from  Equation  26  as 

k  =  |  =  =  k  (0.  947  +  0.  474)  (28) 

s  Z  Z  cp 

or  in  terms  of  the  parameter  G, 


k 

s 


1.0  f  0.  5] 
Z  / 


(29) 


Based  on  these  relationships  between  k  ,  G,  and  kg,  and  noting  that  the 
tire  contact  area.  A,  is  given  by  (13) 

A  =  fDs  (30) 


where 


f  =  coefficient  related  to  the  tire  type  and  tire  deflection 
D=  tire  diameter,  in. 

Bekker's  sinkage  equation  can  be  put  in  the  form 

Z  _  0.  947P 
D  =  GfD3 


(31) 


where 


Z/D  =  3inkage  ratio 


Combining  Equations  24  and  29  leads  to  Boeing's  sinkage  equation  in  the 
form 

Z  _  P  -  0.  5  VGfP8 
D  =  YGfD3 

The  appropriate  Y  values  were  developed  by  Boeing  (5).  By  combining 
Equations  18  and  20,  WES's  sinkage  equation  becomes 


(32) 


~  -  0.003  + 


G(bD)3/* 

P 


6  1 

t 

h 

tJ 

. 

1.  B 


(33) 


In  order  to  compare  these  sinkage  equations,  the  trend  of  the  sinkage 
ratio  versus  increasing  load  was  examined  for  a  20.00-20  type  III  tire  at 
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33-1/3%  tire  deflection  (d).  The  results  of  this  comparison  are  given  in 
Figure  9  for  the  common  soil  parameter  ,  G,  taking  on  values  of  5,  10,  and 
20.  The  G  values  of  5  and  20  could  approximate  a  frictional  soil  in  a  loose 
and  medium  dense  condition  respectively.  The  corresponding  values  of  k 

s 

for  the  Boeing  sinkage  equation  and  k  for  the  Bekker  sinkage  equation  for 

cp 

each  selected  G  value  as  developed  from  Equations  27  and  28  are  also  indi¬ 
cated  in  Figure  9.  This  comparison  is  not  an  attempt  to  study  the  exactness 
of  the  magnitude  of  sinkage,  but  rather  to  examine  trends  with  increasing 
loads.  It  is  evident  however  that  even  though  each  sinkage  equation  utilizes 
the  same  G  value  based  on  the  above  correlation,  that  the  sinkage  ratio  varies 
by  up  to  100%  or  more  between  these  sinkage  equations  under  like  conditions. 
Also  it  is  noted  that  the  trend  of  the  sinkage  ratio  prediction  as  given  by 
Boeing  is  such  as  to  intersect  the  trend  of  the  sinkage  ratio  prediction  as 
given  by  the  Bekker  and  WES  equation. 

One  way  in  which  increased  flotation  (reduced  sinkage)  can  be 
developed  for  pneumatic  tires  is  by  increasing  the  tire  deflection  (decreas¬ 
ing  tire  inflation  pressure).  In  order  to  examine  the  trend  of  this  increased 
flotation  with  increasing  tire  deflection.  Equations  31,  32,  and  33  were 
evaluated  for  a  20.00-20  tire  subjected  to  a  6000  lbt  wheel  load  on  a  G  =  5 
sand  soil.  The  results  of  this  evaluation  are  shown  in  Figure  10.  While  no 
system  has  been  developed  for  rating  increased  flotation  as  related  to  tire 
and  soil  variables,  it  is  evident  from  reference  to  Figures  10  and  5  that  by 
increasing  the  tire  deflection  from  15%  to  30%,  a  decrease  in  the  magnitude 
of  drag  forces  of  between  40%  and  50%  will  occur  for  the  20.  00-20  tire  under 
these  conditions. 

Increased  flotation  can  also  be  developed  by  increasing  the  width  of 
tires.  Figure  11  indicates  the  variation  in  sinkage  ratio  with  increasing 
tire  width  for  a  55  inch  diameter  tire  at  30%  tire  deflection  on  a  sand  soil 
(G  =  5)  as  determined  from  Equations  31,  32,  and  33.  While  the  magnitude 
of  the  sinkage  ratio  is  quite  different  as  determined  from  each  sinkage 
equation,  reference  to  Figures  11  and  5  indicates  that  the  effect  of  varying 
the  tire  -width  from  15"  to  25”  is  to  decrease  the  magnitude  of  drag  by 
approximately  50%  under  these  conditions. 

b)  Sinkage  in  Clay 

Correlation  of  the  soil  parameters  related  to  clay  type  soils  is 
not  as  evident  as  in  the  case  of  sands  and  must  rely  to  a  limited  extent  on 
the  empirical  relationships  established  by  experimental  test  data  as 
related  to  the  CBR  and  cone  index. 

For  clay  type  soils,  experimental  data  indicates  a  range  of  n 
values  in  Bekker’s  equation  of  between  0.  15  to  0.40.  From  Janosi  (12)  a 
relationship  can  be  established  between  the  CBR  at  0.  2  inch  penetration  and 
Bekker's  soil  parameters  k^  and  k  using  n  =  0.  25  as  given  in  Equation  34. 
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Figure  9  Sinkage  Ratio  vs  Wheel  Load,  20.  00-20  Tire 
d  =  33-1/3%  on  Sand  Soil,  G  =  5,  10,  20 


(34) 


CBR=  rS|T+kcpj{0*2)  X450( 

where 

r  =  radius  of  the  CBR  piston. 

For  r  =  1.  97",  Equation  (34)  becomes 

22.  3  CBR  =  1.  02  k  +  k  (35) 

c  cp 

The  results  of  an  extensive  correlation  between  the  CBR  and  the  cone  index 
by  Scala  (14)  indicated  that  for  clay  type  soils,  the  cone  index  at  5  inches 
penetration  can  be  related  to  the  CBR  at  0.  2  inch  penetration  through  the 
approximate  relationship;  Cl  =  50  CBR.  Equation  25  becomes  then 


Cl  =  50  CBR  =  1.625 


(Z+1.5) 


n+1  „n+l 


+  0. 5175  k 
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Z-+2 
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(Z+l.  5)Z 
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co 

n+l| 


(Z+l.  5) 
(n+1)  (n+2) 


(36) 


which  for  n  =  0.  25  and  Z  =  5"  gives 


50  CBR  =3.  77  k  +  1.43  k  (37) 

c  cp 

Equations  35  and  37  can  then  be  used  to  define  the  relationship  between  k  , 

c 

k  ,  and  CBR  as 

cp 

kc  =  7.  9  CBR  (38) 


k  =14.4  CBR 

cp 

Utilizing  the  relationship  of  Scala  (Cl  =  50  CBR)  and  Equation  38,  the  Bekker 
and  WES  sinkage  prediction  equations  can  be  put  in  the  form  of  sinkage  ratio 
equations  as  related  to  the  common  soil  parameter,  CBR.  Bekker's  equa¬ 
tion  becomes 

Z _ P4b4 _ 

D  ~f4D9(CBR)4(7.  9  +  14.4b}4  (39) 

The  sinkage  prediction  equation  of  WES  can  be  written  as 


1  =  0. 
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Equations  39  and  40  were  evaluated  for  a  20.  00-20  type  III  tire  at 
33-1/3%  tire  deflection  for  CBR  soil  values  of  2  and  5.  The  results  of  this 
analysis  are  given  in  Figure  12.  As  for  sands,  there  is  considerable  differ¬ 
ence  in  the  magnitude  of  the  sinkage  ratios  in  the  clay  soil  between  each  of 
the  sinkage  prediction  equations. 
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The  effect  of  increasing  tire  deflections  on  the  flotation  in  clay  soils 
is  shown  in  Figure  13  for  the  same  20.  00-20  tire  for  a  CBR  of  2  soil  and  a 
wheel  load  of  16,  000  lbs.  Reference  to  Figures  5  and  13  show  that  by  in¬ 
creasing  the  tire  deflection  from  20%  to  30%,  an  approximate  50%  to  60% 
decrease  in  the  magnitude  of  drag  forces  will  occur  under  these  particular 
conditions. 

4.  Summary 

"While  the  above  detailed  sinkage  relationships  are  not  a  complete 
survey  of  proposed  load-sinkage  models,  they  do  in  general  typify  past  and 
current  efforts  in  the  development  of  load-sinkage  models  for  soil.  At  the 
present  time,  as  indicated  in  Figures  9  through  13,  and  as  observed  from  a 
review  of  the  mobility  and  flotation  literature,  the  reliability  of  sinkage 
prediction  using  current  load-sinkage  equations  would  appear  to  be  in  the 
range  of  ±50%  to  ±100%.  Since,  as  indicated  in  Figure  5,  the  magnitude  of 
drag  varies  linearly  with  sinkage  {for  other  factors  being  constant),  similar 
uncertainties  would  be  associated  with  estimating  drag  magnitudes. 

In  analyzing  soil  behavior  under  wheel  loads  in  an  effort  to  develop 
an  analytical  expression  for  load-sinkage  based  on  measured  soil  properties 
for  flotation  analysis,  several  factors  are  fundamental  to  the  problem. 

These  factors  include: 

a)  The  analytical  expression  (model)  should  be  relatively  simple 

to  facilitate  its  use  in  predicting  sinkages  as  related  to  soil  type. 

b)  The  resultant  model  should  adequately  describe  the  known  response 
of  soil  and  consequently  to  varying  extents  must  include:  elastic, 
and  plastic  behavior,  and  dynamic  effects. 

c)  Unloading  characteristics  must  be  included  to  analyze  subsequent 
wheel  passage  effects. 

d)  The  resultant  model  should  include  soil  parameters  which  are 
amenable  to  the  existing  soil  experimental  property  measurements. 

While  the  above  summarized  sinkage  equations  satisfy  part  of  these 
factors,  they  are  in  general  tied  to  an  empirical  approach  including  only  the 
surface  variables,  and  are  lacking  in  the  inclusion  of  plastic,  dynamic,  and 
subsequent  loading  effects.  Recent  studies  (15,  16)  have  shown  that  plastic 
behavior  of  soil  is  important  in  analyzing  the  response  of  soil  to  loading. 
Plastic  response  is  a  major  factor  for  contact  stress  levels  normally  en¬ 
countered  from  pneumatic  aircraft  tires. 

Until  sinkage  analysis  includes  the  distribution  of  stress  within  the 
soil  mass  and  incorporates  plastic  and  dynamic  effects  as  related  to  the 
stress  distribution,  increased  reliability  of  sinkage  estimates  will  be  made 
only  through  extensive  full  scale  experimental  testing  over  a  wide  range  of 
variables.  Recent  work  by  Yong  (17)  has  given  some  insight  into  the 
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dynamics  of  soil  response  under  moving  wheel  loads.  Additionally  the  use 
of  finite  element  techniques  as  applied  to  the  contact  element-soil  mass 
problem,  wherein  the  soil  response  can  be  treated  as  an  elasto-plastic 
material,  should  provide  reliable  estimates  of  sinkage  as  related  to  the 
landing  gear,  load,  and  soil  variables. 


SECTION  V 


FLOTATION  PARAMETER 

1.  General 

One  of  the  most  critical  needs  of  aircraft  landing  gear  design 
engineers  is  the  development  of  a  suitable  criteria  for  comparing  the  anti¬ 
cipated  performance  (flotation)  of  landing  gear  systems  on  semi-  and 
unprepared  runways  as  an  aid  in  selecting  the  higher  flotation  system. 

The  results  of  the  drag  and  sinkage  analysis  indicated  that  sinkage  is  an 
extremely  important  parameter  in  analyzing  performance.  A  suitable 
criteria  then  in  developing  a  flotation  parameter  (or  relative  merit  index) 
would  be  to  relate  performance  to  sinkage  and  to  base  comparative  tire 
studies  and  landing  gear  system  studies  on  sinkage. 

The  sinkage  of  an  aircraft  tire  is  related  to  the  load,  tire  charac¬ 
teristics,  and  soil  conditions.  Since  only  a  comparison  of  different  tires 
performance  is  to  be  made,  tho  soil  conditions  will  remain  constant  for 
comparison  purposes  which  simplifies  the  problem.  Ideally,  the  deter¬ 
mination  of  the  sinkage  would  be  based  on  total  sinkage  as  developed  from 
an  elasto-plastic -dynamic  theory  which  adequately  reflects  the  deform- 
ability  characteristics  of  soil.  As  indicated  in  Section  IV,  such  a  solution 
is  not  available  at  the  present  time.  In  order  to  develop  the  concept  ~f  a 
flotation  parameter  based  on  sinkage  and  to  indicate  how  comparative  tire 
studies  could  be  made  using  nomographic  charts,  the  existing  theories  of 
an  elastic  half  space  subjected  to  dynamic  loads  were  extended  for  applica¬ 
tion  to  aircraft  operations  on  unprepared  runways.  The  development  of  the 
flotation  parameter  then  incorporates  the  numerous  tire  variables  and  is 
based  on  the  criteria  of  sinkage  where  sinkage  is  determined  based  on 
elastic  theory. 

As  techniques  become  available  which  permit  reliable  predictions 
of  total  sinkage  (elastic  plus  plastic)  the  concept  of  the  flotation  parameter 
which  utilizes  this  total  sinkage  can  become  an  effective  design  criteria  for 
landing  gear  systems. 

2.  Development  of  the  Flotation  Parameter 

a)  Tire  Characteristics 

The  numerous  tire  variables  which  influence  sinkage  have  been 
previously  summarized  in  Figure  2.  Analysis  of  existing  experimental 
data  (1)  has  indicated  that  ply  rating  and  tread  pattern  do  not  significantly 
influence  tire  performance.  The  tire  inflation  pressure,  initial  tire  stiff¬ 
ness,  and  wheel  load  all  influence  the  tire  contact  area,  contact  geometry, 
and  contact  stress  distribution.  If  tire  deflection  (in  per  cent)  is  analyzed 


in  relation  to  contact  area  then  the  complex  interaction  between  tire  stiff¬ 
ness,  inflation  pressure,  and  load  can  be  simplified.  Figure  14  shows  the 
variation  of  normalized  contact  area  (Aj/Ag)  with  normalized  tire  deflection 
(di/d2).  In  preparing  Figure  14,  data  giving  tire  contact  area  at  several 
deflections  for  a  particular  tire  was  used.  Contact  area  and  tire  deflections 
were  normalized  so  that  data  from  several  sources  (1,  7,  18)  could  be 
analyzed.  Tire  diameters  between  30"  and  70"  are  included  which  is  the 
usual  range  for  aircraft  tires.  Reference  to  Figure  14  indicates  that  an 
essentially  linear  relationship  exists  between  normalized  contact  area  and 
normalized  tire  deflection  for  the  range  of  normalized  tire  deflections 
considered  (0.  3  to  2.  5).  This  simplified  relationship  will  permit  comparison 
of  tire  performance  at  different  deflection  levels. 

Information  is  notably  lacking  on  the  distribution  of  contact  stress 
between  tires  and  soil  at  different  inflation  pressures.  Analysis  of  the 
limited  information  available  (18)  indicates  that  at  lower  inflation  pressures 
(high  tire  deflection),  the  contact  stress  distribution  tends  to  be  uniform 
while  at  high  inflation  pressures,  the  contact  stress  distribution  is  parabolic 
with  the  peak  stress  occurring  immediately  beneath  the  center  line  of  the 
tire.  A  region  of  transition  for  contact  stress  distribution  will  naturally  lie 
between  these  two  extremes.  Existing  data  does  not,  however,  permit  a 
quantitative  look  at  this  transition  region.  The  flotation  parameter  nomo¬ 
graph  was  developed  for  the  single  case  of  a  uniformly  distributed  load.  As 
additional  information  becomes  available  on  contact  stress  distribution  vs. 
tire  deflection  in  the  transition  region,  a  correction  factor  can  be  incorpora¬ 
ted  into  the  flotation  parameter. 

The  contact  geometry  between  tires  and  soil  has  not  been  determined 
in  detail.  Table  III  summarizes  a  limited  amount  of  data  from  two  sources 
(8,  18),  which  indicates  the  contact  geometry  to  be  elliptical  in  shape, 

TABLE  in 


Elliptical  Contact  Geometry 


Inflation 

Major 

Minor 

Ratio 

Tire 

Pressure  (psi) 

Axis=a 

Axis=b 

=  a/b 

Soil  Type 

11.00-20, 

12  PR 

15 

18.  5"- 

12.  5" 

1.48 

Sand 

11.  00-20, 

12  PR 

15 

19.  0" 

12.  5" 

1.  52 

Clay 

11.  00-20, 

12  PR 

30 

12.  0" 

7.  0" 

1.72 

Rigid 

vO 

o 

o 

* 

h- * 

8  PR 

15 

9.  3" 

6.  3" 

1.46 

Rigid 

As  indicated  in  Table  III,  the  contact  shape  tended  to  be  elliptical 
with  ratios  of  the  major  to  minor  principal  axes  ranging  from  1.  4  to  1.  7. 

This  limited  data  indicates  that  the  contact  geometry  does  not  differ  signi¬ 
ficantly  from  a  circular  contact  area,  and  consequently,  since  only  maximum 
vertical  sinkage  is  required,  circularly  loaded  areas  were  used  in  the 
development  of  the  flotation  parameter. 
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REFERENCE  TIRES 


b)  Ground  Loading 

The  determination  of  the  load  time  history  at  a  point  in  a  semi¬ 
infinite  half  space  (soil  medium)  due  to  aircraft  loading  is  of  considerable 
importance  in  analyzing  the  load-sinkage  relationship  for  moving  aircraft 
wheels.  The  dynamic  ground  sinkage  depends  not  only  on  the  peak  load  but 
also  on  the  shape  and  duration  of  the  load  pulse.  The  ground  surface  load¬ 
ing  can  be  thought  of  as  having  two  components,  vertical  and  horizontal 
(shear  or  drag)  which  influence  the  sinkage  of  the  aircraft  tire  into  the  soil 
(see  Figure  3).  Table  IV  summarizes  the  relative  magnitude  of  each  of  these 
components  in  defining  sinkage  for  the  different  types  of  aircraft  loadings 
on  soil. 

TABLE  IV 

Influence  of  Component  Loading  on  Sinkage 


Components  of  Resultant  Ground  Loading 


Type  of  Loading 

Vertical 

Horizontal 

Landing  (Impact) 

Major 

Minor 

Taxiing  (Dynamic) 

Region  2 

M  a  j  or 

Minor 

Takeoff  (Dynamic) 

Velocity  Range 

Major 

Minor 

Braking  (Dynamic) 

Major 

Major 

Reference  to  Figure  5  indicates  that  in  most  instances  the  ratio  of  the  drag 
force  to  the  vertical  force  is  one-quarter  or  less  for  taxiing  and  takeoff 
operations.  Since  the  flotation  parameter  is  to  be  utilized  for  comparative 
purposes  only,  and  noting  that  the  major  factor  in  vertical  sinkage  is  the 
vertical  load,  the  flotation  parameter  was  developed  based  on  vertical  load¬ 
ings.  Little  information  is  available  concerning  the  shape,  duration,  and 
magnitude  of  the  ground  loading  for  landing  situations  and  consequently  the 
flotation  parameter  nomograph  is  related  only  to  taxiing  and  takeoff  operations. 
The  impact  landing  and  braking  loading  conditions  are  to  be  studied  in  a  sub¬ 
sequent  phase  of  the  research  program. 

Considerable  information  is  available  on  the  load-time  histories 
at  a  point  in  a  semi-infinite  half  space  under  moving  wheel  loads  from 
research  work  in  flexible  pavements  and  soil  mechanics.  Experimental 
measurements  of  load-time  histories  by  use  of  load  cells  embedded  in 
flexible  pavements  and  in  soil  for  rolling  wheels  provides  an  excellent 
source  for  defining  the  expected  ground  loading  conditions  in  the  path  of 
an  aircraft  wheel  for  taxiing  and  takeoff. 

Figure  15  shows  several  experimentally  measured  load-time 
histories  (19,  20)  which  are  typical  of  embedded  load  cell  measurements. 

For  comparison  purpose,  a  half  si?ie  wave  has  been  superimposed  on  the 
actual  load-time  trace.  There  is  considerable  similarity  between  the  half 
sine  wave  and  the  actual  load-time  trace.  Based  on  such  comparisons,  it 
has  been  common  practice  to  assume  that  the  load-time  history  at  a  point 
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in  a  soil  medium  or  flexible  pavement  under  moving  wheel  loads  can  be 
adequately  approximated  by  a  half  sine  wave  (compressional  loading).  It  has 
further  been  shown  that  the  di  ration  of  the  half  sine  wave  loading  is  directly 
related  to  the  speed  of  the  traversing  wheel.  Table  V  summarizes  wheel 
velocity -loading  duration  information  currently  available  as  developed  from 
research  work  in  soil  mechanics  and  flexible  pavements. 


Horizontal 
Wheel  Velocity 
(knots) 

60 

30 

4 


TABLE  V 

Velocity-Duration  Relationship 

Duration  of 
Compressional  Wave 
(sec. ) 

0.  02  to  0.  05 
0.  11  to  0.  14 
0.  85  to  1.  00 


Reference 


21,  22 
22,  23 
22,  23 


Based  on  the  above  detailed  experimental  information,  the  load¬ 
time  history  at  a  point  in  the  soil  medium  subjected  to  dynamic  loading  by 
an  aircraft  wheel  will  be  represented  by  a  half  sine  wave  loading  pulse.  On 
a  preliminary  basis,  Table  V  will  be  utilized  for  relating  horizontal  ground 
velocity  to  the  loading  duration. 


c)  Dynamic  Soil  Response 


The  response  of  a  semi-infinite  soil  mass  under  dynamic  loading 
has  been  determined  on  the  assumption  of  the  soil  being  a  homogenous 
elastic  medium  (24,  25).  Recent  theoretical  work  by  Ri chart  and  Whitman 
(26,  27)  has  demonstrated  that  for  determining  the  relationship  between  load 
and  soil  response  (vertical  sinkage),  that  the  dynamic  behavior  of  the  system 
can  be  adequately  represented  as  a  single  degree  of  freedom  lumped  parameter 
system  as  indicated  in  Figure  16.  The  viscous  damping  component  as 
represented  by  the  dash  pot  in  the  model  primarily  reflects  the  damping 
action  due  to  the  dispersion  of  the  elastic  waves  throughout  the  soil  medium. 
The  differential  equation  of  motion  for  the  model  shown  in  Figure  16  for 
uniform  loading  is  given  by 

mZ  +  CjcZ  +  kxkZ  =  P(t)  (41) 


where 


m 

c 


k 


equivalent  mass  of  gear-wheel  system 

kr 
_ o 

V 

s 

irGr 
_ _ o 

(1-H) 
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and 


=  Poisson's  ratio 

r  =  radius  of  circularly  loaded  area 
o 

G  =  soil  shear  modulus 
Cijkj  =  parameters  defined  in  Appendix  II 


P(t  ) 


Figure  16  Half  Space  Mathematical  Model  (Elastic  System) 


The  forcing  function  P(t)  is  the  time  dependent  vertical  force.  The  solution 
of  Equation  (41)  for  a  steady  state  harmonic  loading, 


where 

is  given  fcy 
where 


P(t)  =  P  sinujt 
o 


Pq  =  peak  amplitude  of  harmonic  force 
0)  =  angular  frequency  of  steady  state  motion 
P„ 


(42) 


Z(t)  =  —  M  sin  (wt  +  cp) 
1 


(43) 


M  = 


V?i-Bao3)3  +  (cia0)! 


_  l-U  m 

B  =  mass  ratio  =  •  -'r~ ■  y 

o 
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and 


p  =  mass  density  of  soil 


a  =  frequency  ratio  =  uur  ^p/G* = 


(Dr 
■»  o 


V 


/ci a  * 

„  _  t,_~i  (  °) 

<P"tan  1-Ba  3 
o 


V  =  shear  wave  velocity 
s 


As  indicated  previously,  the  aircraft  tire  is  moving  with  a  hori¬ 
zontal  velocity  across  the  soil  medium.  Consequently,  the  loading  function 
P(t)  is  a  transient  pulse  in  the  form  of  a  half  sine  wave  of  duration,  t^(see 

Figure  15).  Presently  available  solutions  are  for  steady  state  harmonic 
loading  functions  as  indicated  by  Equation  43;  however,  Richart  (26)  has 
suggested  a  method  of  solution  for  transient  loadings  by  use  of  Fourier 
analysis.  By  determining  a  Fourier  series  of  loading  functions  which 
approximate  the  transient  loading,  and  using  the  principle  of  superposition 
to  add  together  the  individual  displacement  solutions,  the  sinkage  response 
is  given  in  the  form 

N 

P  <-— > 

Z(t)  =“  L  M  sM«)t  +  o>  -  Y  )  (44) 

k  1 — ■  n  n  n  n 

n=0 


where 


Y  =  phase  relationship  between  the  individual  frequency 
components 


Computer  programs  were  written  to  determine  the  maximum 
sinkage  under  a  half  sine  wave  transient  loading  for  the  variables  of  load 
duration  (t^),  soil  properties  (|i,  G,  p),  mass  (m),  and  peak  amplitude  of 

load  (Pq).  The  results  were  then  used  to  develop  the  flotation  parameter 

nomograph.  The  full  details  of  the  transient  loading  solution,  together  with 
the  computer  program,  are  given  in  Appendix  II. 


3.  Flotation  Parameter 

Flotation  is  defined  as  the  ability  of  the  soil  to  support  the  load 
induced  through  the  contacting  element.  The  flotation  parameter  (F)  is 
designed  to  absorb  the  numerous  tire  variables  into  a  single  parameter  for 
comparative  tire  (or  contacting  element)  studies  based  on  flotation.  The 
criteria  for  comparison  of  flotation  in  determining  flotation  parameters  is 
sinkage.  The  flotation  parameter  nomograph  shown  in  Figure  17  is  based 
on: 

(1)  The  complex  interrelationship  of  the  tire  variables  can  be 
simplified  through  the  relationship  of  Figure  14. 
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(2)  The  contact  stress  is  assumed  to  be  uniform.  * 

(3)  The  contact  area  is  assumed  to  be  circular.  * 

(4)  The  loading  at  a  point  in  the  soil  mass  due  to  the  moving  wheel 

is  in  the  form  of  a  half  sine  wave  transient  load  of  duration,  t ,. 

d 

(5)  The  soil  mass  is  assumed  to  be  elastic. 

(6)  The  mass  of  the  gear-wheel  is  incorporated  through  the  mass 
ratio  parameter,  B.  If  B  =  0,  then  only  a  forcing  function  exists 
at  the  ground  surface. 


Example  of  Use  of  Nomograph 

Problem  1  Compare  the  relative  flotation  characteristics  of  the 
56  x  16  tire  to  the  40  x  14  tire  for  both  tires  operating  at  30%  tire  deflection 

(d).  , 

Known  Information 

i 

Clay  soil,  shear  wave  velocity,  V  =  250  fps, 

horizontal  wheel  velocity  =  25  knots,  f,  =  0.  20  sec, 

d 

the  56  x  16  tire  has  a  known  bearing  area  of  340  in.  8  at  d  =  40% 
the  40  x  14  tire  has  a  known  bearing  area  of  97  in.  8  at  d  =  20% 

Use  B  =  0 

Solution  by  Nomograph 

56  x  16  tire,  D  =  55.  7" 

(1)  dj/dg  =  40/30  =  1.  33 

(2)  V  •  t ,  =  250  x  12  x  0.  2  =  600 

s  d 

(3)  Intersect  340  in.  8  (Ax)  with  dT  /d2  =  1.  33  which  gives 

r  at  d  =  30%  of  8.  8". 
o 

r 

(4)  Intersect  r  =  8.  8"  with  V  •  t ,  of  600  to  get  — —  =  0.  015 

o  s  a  V  •  t 

8  d 

(5)  Move  horizontally  to  intersect  B  =  0  curve  which  gives 
M  *  1.0 

(6)  Intersect  M  =  1.  0  with  rQ  =  8.  8"  to  get  the  flotation 

parameter  F_,  =2.5 

56  x  16 


Adjustments  can  be  made  to  reflect  more  exact  conditions,  if  experimental 
data  becomes  available  which  indicates  such  requirements. 


4 

I 
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40  x  14  tire,  D  =  39.  3" 

(1)  di/da  =||=  0.67 

(2)  V  *  t  =  250  x  12  x  0.  2  =  600 

s  d. 

(3)  Intersect  97  in.  8  (Aj)  with  d^dg  =  0.  67  which  gives 

r  =  7.  7" 
o 

(4)  Intersect  r  =7.  7"  with  V  't,  of  600  to  get 

o  s  d  & 

v-ftT  =  0-013 

s  d 

(5)  Move  horizontally  to  intersect  B  =  0  curve  which  gives 
M  =  1.0 

(6)  Intersect  M  =  1.  0  with  r  =7.  7"  to  get  the  flotation 

o 

parameter  F  .  ,  =  2.  9 

r  40  x  14 

The  56  x  16  tire  has  higher  flotation  characteristics  since  it  has  less 
sinkage, 

F56  x  16  <  F40  x  14 

Relating  the  flotation  to  drag  by  use  of  Figure  5,  would  indicate  an 

even  higher  relative  flotation  rating  for  the  56  x  16  tire  since  drag  varies 

with  the  sinkage  ratio,  (Z/D),  and  D..  ,,  is  larger  than  D.„  ,  thus: 

56  x  16  40  x  14 

F  j  ,  *>.*-«-  ?§.  <  %  *  J4,  g  F  1 

56x16  D  ,  D .  ..  ’40x14 

56  x  16  40  x  14 


where 


F.  I.  =  flotation  index 


4.  Summary 

The  nomograph  shown  in  Figure  17  is  an  attempt  to  develop  the  con¬ 
cept  of  a  flotation  parameter  in  order  to  permit  comparative  flotation  studies 
among  landing  gear  systems.  By  using  Figure  5  with  the  flotation  parameter 
nomograph,  comparison  studies  can  also  be  made  based  on  the  flotation 
index.  It  does  not  provide  completely  adequate  comparisons  at  the  present 
time,  since  the  sinkage  is  based  solely  on  elastic  theory  while  total  sinkage 
of  an  aircraft  tire  involves  both  elastic  and  plastic  deformations.  As 
reliable  techniques  become  available  for  predicting  sinkage,  the  nomograph 
can  be  suitably  modified  to  reflect  these  techniques  and  can  provide  the 
landing  gear  design  engineer  with  considerable  flexibility  in  analyzing  flota¬ 
tion  with  the  following  objectives: 


44 


(1)  Comparative  tire  flotation  studies  where  the  soil  and  loading 
conditions  are  constant  (see  example  problem  #1). 

(2)  Comparison  of  performance  of  a  tire  on  different  soil  types. 

(3)  Design  of  landing  systems  and  contacting  elements  for  reduced 
sinkage. 

(4)  Selection  of  tires  to  reduce  sinkage. 

(5)  Determination  of  increased  flotation  with  increasing  tire 
deflections. 

(6)  Comparative  studies  of  drag  magnitudes  for  the  region  two 
velocity  range  (see  example  problem  #1). 

(7)  Development  of  operational  criteria  for  aircraft  at  specific 
landing  sites. 

These  comparative  studies  conducted  by  use  of  a  flotation  parameter 
nomograph  based  on  total  sinkage  will  lead  ultimately  to  a  comparison  of  the 
relative  merits  of  landing  gear  configurations  with  regard  to  their  flotation 
capability  on  semi-  and  unprepared  soil  (relative  design  merit  index). 


SECTION  VI 


AIRCRAFT  OPERATION  CAPABILITY  (SUITABILITY  OF  SITE) 

1,  Failure  Criteria 

Determination  of  the  capability  of  aircraft  to  operate  on  semi-  and 
unprepared  sites  is  generally  made  based  on  performance  and  is  expressed 
as  the  number  of  coverages  that  can  be  carried  out  previous  to  failure 
(coverage  is  defined  as  sufficient  passes  of  load  tires  in  adjacent  tire  paths 
to  cover  a  given  width  of  surface  area  one  time).  Current  Air  Force  stan¬ 
dards  specify  the  failure  criteria  as  occurring  when  permanent  rutting  of 
clay  soil  exceeds  3  inches.  The. Army  (9),  in  a  recent  experimental  research 
investigation  of  aircraft  tires  operating  on  unprepared  surfaces,  utilized  a 
failure  criteria  of  3  inches  permanent  deformation  or  1.  5  inches  elastic 
deformation.  The  results  of  this  research  led  to  a  modified  "Unsurfaced 
Soil  Strength  Requirements"  nomograph  for  determining  the  allowable 
number  of  aircraft  tire  coverages  on  forward  airfields.  For  reference 
purposes  this  nomograph,  which  includes  the  variables  of  load,  tire  infla¬ 
tion  pressure,  and  CBR,  is  shown  in  Figure  18.  For  multiple  wheel  arrange¬ 
ments  an  equivalent  single  wheel  load  is  determined,  based  on  Figure  19.  to 
account  for  the  effects  of  adjacent  wheels  as  related  to  tire  spacing. 

Several  limitations  exist  in  the  use  of  the  strength  requirements 
nomograph.  The  relationships  developed  are  based  on  experimental  data 
from  tire  tests  primarily  on  "buckshot"  clay  soil.  The  recent  modifica¬ 
tions  were  made  from  test  data  on  Buckshot  clay  which  exhibits  a  high  degree 
of  elastic  response.  Consequently,  as  noted  by  the  Air  Force  (28),  the  use 
of  the  nomograph  for  determining  the  number  of  operations  is  applicable  only 
to  heavy  clay  soils.  Also,  the  distribution  of  loads  used  in  defining  coverage 
by  the  Army  (9)  often  do  not  apply  in  forward  airfield  situations  since  the 
limited  width  airfields  restrict  aircraft  operations  to  the  centerline  of  the 
airstrip.  Recent  returning  pilots  from  Vietnam  who  have  operated  C-123 
and  C-I30  aircraft  on  unprepared  airstrips  report  satisfactory  operations 
on  unprepared  airstrips  with  rut  depths  up  to  five  inches.  Additionally, 
pilots  have  indicated  the  tendency  to  continually  land  aircraft  in  the  center- 
line  rut  due  to  the  limited  airfield  widths  normally  encountered  in  forward 
areas.  The  present  criteria  also  does  not  include  braking,  turning,  or 
landing  impact  effects. 

2.  Unified  Approach 

The  problem  of  determining  aircraft  capability  to  operate  at  a  site 
for  aircraft  operations  clearly  involves  the  relationship  between  permanent 
rut  depth,  roughness,  and  steering  ability  in  addition  to  required  runway 
lengths.  The  combination  of  initial  roughness  with  subsequent  permanent 
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INCREASE  IN  SINGLE-WHEEL  LOAD  FOR  EACH 
ADJACENT  WHEEL,  PERCENT* 
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C-C  TIRE  SPACING,  RADII 


*  INCREASE  IN  LOAD  ON  A  SINGLE  WHEEL  OF  A 
MULTIPLE-WHEEL  GEAR  TO  ACCOUNT  FOR 
EFFECTS  OF  ADJACENT  WHEELS  OF  THE 
MULTIPLE-WHEEL  GEAR  IN  ARRIVING  AT  AN 
EQUIVALENT  SINGLE-WHEEL  LOAD. 


Figure  19  Equivalent  Single-Wheel  Load  Adjustment 
Curve  for  Unsurfaced  Soils 
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deformations  defines  the  roughness  of  the  airstrip  for  subseqiient  operations. 
Figure  20  shows  the  interrelationship  of  these  factors  in  th3  development  of 
a  unified  approach  to  aircraft  operations  capability  and  flotation.  Since  dif¬ 
ferent  types  of  aircraft  undergo  different  dynamic  structural  responses  when 
traversing  the  same  rough  terrain,  the  suitability  should  be  related  to 
specific  types  of  aircraft.  Until  either  additional  full  scale  testing  results 
become  available  on  accepted  tolerances  of  aircraft  to  rough  terrain  or  until 
more  accurate  sinkage  prediction  techniques  are  developed,  the  capability  of 
an  aircraft  to  operate  at  a  site  cannot  be  fully  defined. 
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Figure  20  Unified  Approach:  Flotation  Parameter- 
Aircraft  Operation  Capability 
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SECTION  VII 


CONCLUSIONS  AND  RECOMMENDATIONS  FOR  RESEARCH 
1.  Conclusions 

In  general,  the  results  of  the  initial  phase  of  the  research  effort 
including  the  literature  review,  have  shown  that  the  effects  of  such  variables 
as  soil  type,  speed,  tire  size,  and  other  factors  have  not  been  completely 
defined  in  relation  to  the  flotation  of  aircraft  tires  on  soil.  Those  variables 
which  significantly  effect  the  landing  gear-soil  interaction  phenomena  have 
been  identified.  A  review  of  existing  flotation  and  mobility  data  has  shown 
that: 

(1)  There  are  at  least  three  velocity  regions  for  which  the  drag 
ratio  exhibits  distinct  characteristics. 

(2)  In  the  region  two  velocity  range,  the  drag  ratio  (R/P)  is  directly 
proportional  to  the  sinkage  ratio  (Z/D)  for  single  tires. 

(3)  Information  is  lacking  on  those  variables  which  define  the  drag 
ratio  response  La  the  region  three  velocity  range. 

(4)  Little  information  is  available  which  shows  the  effect  of  landing 
impact,  braking  (including  antiskid  system  operation),  and 
turning,  on  flotation  and  operation  capability. 

Since  sinkage  has  been  shown  to  be  a  significant  variable  in  defining 
flotation,  a  review  of  the  existing  sinkage  studies  indicated  three  currently 
used  sinkage  prediction  methods  (models).  The  results  of  the  comparative 
study  between  these  sinkage  prediction  equations  under  identical  soil, 
geometry,  and  loading  conditions  showed  considerable  difference  in  the 
magnitude  of  sinkage  between  methods.  While  the  comparative  study  was 
based  on  a  correlation  of  soil  parameters  (kc,  k  ,  n,  C.I.  ,  G,  and  CBR), 

the  result  of  the  comparison  together  with  the  review  of  existing  literature 
indicates  that  the  accuracy  of  sinkage  prediction  utilizing  present  techniques 
is  in  the  range  of  ±50%  to  ±100%.  In  spite  of  the  marked  differences  in  the 
magnitude  of  sinkages  predicted  by  the  different  equations  for  similar 
conditions,  all  the  sinkage  analysis  equations  indicate  clearly  the  quantitative 
changes  in  flotation  (sinkage  and  drag)  with  changes  in  the  variables  of  tire 
deflection  (d),  tire  diameter  (D),  and  tire  widrh  (b).  For  example,  by 
increasing  the  tire  deflection  of  a  20.  00-20  tire  from  20%  to  30%  in  a  clay 
soil  (CBR=2),  an  approximate  decrease  of  50%  results  in  the  magnitude  of 
drag.  Until  more  accurate  sinkage  prediction  techniques  are  developed, 
particularly  in  defining  the  deformability  response  of  soil  to  dynamic  loads, 
suitable  prediction  of  flotation  response  for  determining  the  capability  of 
aircraft  to  operate  from  a  site  cannot  be  made. 
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In  order  to  meet  the  requirements  of  landing  gear  design  engineers 
for  the  development  of  high  flotation  tires  and  landing  gear  systems  based  on 
a  comparative  evaluation,  the  concept  of  a  flotation  parameter  and  flotation 
index  which  ultimately  should  lead  to  a  relative  design  merit  index  was 
developed.  While  the  present  flotation  parameter  utilizes  elastic  theory  in 
defining  soil  response,  as  sinkage  prediction  techniques  become  available 
which  include  elastic  and  plastic  deformation  (total  sinkage),  the  flotation 
parameter  nomograph  can  be  suitably  modified  to  reflect  total  sinkage.  The 
resulting  modified  nomograph  could  be  a  powerful  tool  in  the  comparative 
evaluation  of  landing  gear  systems  for  operation  on  semi-  and  unprepared 
soil. 


2.  Recommendations  for  Research 

The  review  of  the  landing  gear-soil  interaction  problem  has  shown 
several  research  deficient  areas.  These  required  research  areas  are 
summarized  as  follows: 

a;  Development  of  an  accurate  and  reliable  sinkage  model  incorporat¬ 
ing  load,  geometry,  and  soil  conditions.  The  extension  of  present  techniques 
in  finite  element  analysis,  with  a  suitable  yield  criterion  for  soil  to  reflect 
the  plastic  deformation  of  soil,  for  application  to  the  aircraft  tire-soil  inter¬ 
action  problem  should  lead  to  a  more  reliable  sinkage  prediction  technique. 

In  order  to  effectively  develop  the  sinkage  analysis,  suitable  correlations 
should  be  made  between  the  sinkage  prediction  technique  and  experimental 
data  from  a  controlled  sinkage  test  program. 

b.  Further  investigation  of  the  drag  ratio-velocity  spectrum  for  the 
region  three  velocity  range  and  higher  velocities.  Efforts,  both  theoretical 
and  experimental,  should  be  directed  at  the  development  of  a  drag  ratio 
equation  as  related  to  the  significant  variables. 

c.  More  reliable  determination  of  twin  (spacing)  and  tandem  effects 
as  related  to  flotation  for  aircraft  wheels  on  soil  (sand  to  clay). 

d.  Evaluation  of  the  influence  of  braking,  turning,  and  landing 
impact  effects  on  surface  roughness  and  flotation. 

e.  Extension  of  the  present  flotation  parameter  (F)  and  flotation 
index  (F.  I. )  to  a  nomograph  based  on  total- sinkage  for  comparative  studies 
and  rating  of  present  and  proposed  aircraft  tires  and  landing  gear  systems. 

On  a  more  limited  basis,  investigation  into  the  following  areas  could 
provide  considerable  additional  insight  into  the  interaction  response  between 
soil  and  aircraft  tires: 

f.  Tire  contact  pressure  distribution  and  contact  geometry  effects 
on  flotation  in  soil. 

g.  Hydroplaning  phenomena  as  related  to  tire  and  soil  conditions. 

h.  Additional  evaluation  of  the  failure  criteria  for  defining  the 
capability  of  aircraft  to  operate  from  a  site. 
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APPENDIX  I 


DRAG-SINKAGE  EXPERIMENTAL  DATA 
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APPENDIX  II 

TRANSIENT  LOADING  SINKAGE  ANALYSIS, 
COMPUTER  PROGRAM 
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The  steady  state  response  of  a  rigid  body  in  contact  with  a  semi¬ 
infinite  elastic  medium  has  been  analyzed  by  a  number  of  investigators 
(24,  25,  26).  While  there  are  four  modes  of  motion  (vertical,  torsion, 
horizontal,  and  rocking)  for  a  rigid  body  with  a  circular  contact  area,  only 
the  vertical  mode  is  of  interest  at  this  time  to  the  flotation  problem.  For  a 
rigid  body  resting  on  a  semi-infinite  medium  and  using  the  assumptions  of  a 
circular  contact  area,  a  homogeneous,  isotropic,  and  elastic  medium,  and 
onlY  normal  stresses  transferred  at  the  interface,  the  vertical  displacement 
(sinkage)  of  the  rigid  body  with  a  uniform  contact  stress  for  a  harmonic 
loading  of  the  form 

p(t)  =  P  e  (45) 


is  given  by 


where 


rr _ o  .  .  iu)t 

z  gT  {fl  +  lfs)e 

o 


=  dimensionless  functions  dependent  on  Poisson's 

ratio  (p),  and  the  frequency  ratio  (a  ) 

o 

Pq  =  peak  amplitude  of  harmonic  force 

(1)  =  angular  frequency  of  steady  state  motion 
t  =  time 

G  =  soil  shear  modulus 


»-  y/T 


where 


r  =  radius  of  circularly  loaded  area 
U)r 

o 

ao  "  V 

s 


Vg  =  shear  wave  velocity  in  soil 


The  values  of  fj  and  f3  are  given  in  Table  XIII  for  a  uniform  contact  stress. 

More  recently  the  work  by  Lysmer  and  Richart  (26)  has  shown  that 
by  defining  the  quantities 


"Va 

3  o 

''  F,  8  +  F-s 
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where 


Fx  =  irfi/ti-u),  F3  =  nfa/(l-|i) 

An  analogy  can  be  drawn  between  the  response  of  the  half  space  model  and 
a  simple  single  degree  of  freedom  lumped  parameter  system  (see  Figure  16). 
Noting  that  the  differential  equation  of  motion  for  a  weightless  rigid  body  is 
given  by 

1c  X* 

c,  — -  Z+k,kZ  =  P  e1U)t  (49) 


- -  Z  +  k,  kZ  =  P  e 

V  o 

s 


where 


"(1-4) 

and  that  by  adding  the  inertia  term  (m  t  0),  the  differential  equation  of 
motion  becomes 

kr 

ml  +  Cl— -  Z  +  k^Z  =  P  e1Wt  (50 

1  V  1  o 

s 

The  values  of  c1  and  kx  as  developed  from  equations  47  and  48  are  given  in 
Table  XIII.  The  solution  of  Equation  50  for  a  forcing  function 

P(t)  =  P  siniwt  (51 

o 


Z  =  — -  M  sin  (u)t  +  cp) 


where 


'(k-t-Ba^)3  +  (cxa  )* 


B  =  mass  ratio  = 


-  u  m 

7T  or  3 


cp  =  tan' 


-c  a 
i  o 

k,  -  Ba  s 


For  transient  loading  phenomena,  Lysmer  and  Richart  (26)  suggested 
a  Fourier  series  analysis  utilizing  steady  state  solutions.  Figure  21  shows 
a  vertical  loading  in  periodic  form.  The  loading  function  is  a  sequence  of 
identical  pulses  of  alternating  sign  and  duration,  t^.  The  pulses  are  uniformly 
spaced  and  periodic  (2jt^l.  By  selecting  the  integer,  j,  large  (large  spacing 
between  the  pulses),  the  system  will  undergo  free  vibrations  in  response  to 
each  pulse.  These  free  vibrations  will  damp  out  (large  damping  due  to 
energy  dispersion  in  the  half  space)  previous  to  the  loading  by  a  subsequent 
pulse.  In  effect  then  the  system's  response  will  be  as  if  its  loading  were  a 
transient  pulse.  In  this  case  the  transient  pulse  has  been  selected  as  a  half 


'  - 


sine  wave  of  duration,  t^,  and  peak  intensity,  P^.  Putting  the  loading 

function  in  nondimensional  form  as  indicated  in  Figure  22,  the  load  pulse 
can  be  defined  as: 

P(t)  =  P  co(t/t ,)  (56) 

o  a 

Reference  to  Figure  21  indicates  that  the  loading  function  is  of  the  form: 

f(t-  jtd)  =  -f(t)  (57) 

Expanding  the  loading  function  (Equation  56)  in  a  Fourier  series,  letting 

u  =  t/t ,,  and  noting  the  relationship  of  Equation  57  leads  to  the  following 
d 

equations  for  the  Fourier  series  coefficients 


a  =  — 
«  J 


b  =  — 
n  J 


-j  +  1  .1 

/  .  co(u)  cos  du  +  /  cp(u)  cos—-  du 

J~J  J  o  1 

'  -j+1  1 

/  co(u)  sin  ~ —  du  +  /  co(u)  sin— ™  du 

J  o  J 


which  reduces  to 


a  =  —  j  co(u)  cos 

n  j  Jo 

2  1 

b  =  —  /  co(u)  sin 

n  j  n 


(2n+l) 

i  - 

(2n+l)j 


TTU  QU 


TTudu 


(58) 

(59) 

(60) 

(61) 


The  Fourier  series  developed  from  this  loading  function  is  then 

N 

p(t)  =  p  y  c  cos 

o  L-  ^  n 


where 


and 


n=0 


(2n+l)Trt 


i  t. 


-  v 


n 


2  +b5 


c  =  \/a‘ 
n  \/  n  n 

b 

Y  =  tan-1  — 
n  a 

n 


N  is  a  large  integer. 


(62) 


Using  the  principle  of  superposition  as  applied  to  Equation  52  for 
linear  systems  gives  the  following  series  for  the  vertical  displacement 
(sinkage) 
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where 


p  N  * 

Z(t)  =  V  M  cos  irfc  +  cd  -  Y  | 

k  £— •  n  j  t ,  n  n| 

n=0  J  d  J 


M  =, 
n 


(k,  -  Ba  2)2  +  (c  a  ): 

in  on  m  on 


*n  =  tan_1  P^2  <66> 

in  on/ 

The  values  of  c  and  k  can  be  determined  from  Table  XIII  based  on  a 

min  on 

As  indicated  by  Lysmer  and  Richart  (26)  most  of  the  response  is 
developed  from  the  lower  frequency  contents  and  consequently  a  good  approxi¬ 
mation  to  the  sinkage  can  be  obtained  for  reasonably  small  values  of  N.  An 
approximation  to  the  number  of  terms  required  is  given  by  (26) 

N>.f  ‘ 11  V\  («> 

(2irr  J  Be  ) 

o  T 

where  e  defines  the  required  accuracy. 

The  determination  of  the  sinkage,  Z,  through  the  series  as  expressed 
by  Equation  63  requires  the  evaluation  of  a  and  b  by  the  integral  form  of 

Equations  60  and  61.  Since  in  the  general  case  the  load  pulse  as  described 
by  cc(t/td)  is  not  always  a  shape  which  can  be  described  by  a  closed  form 

equation,  the  loading  function  was  specified  as  a  segmented  function  as  shown 

in  Figure  22  and  consequently  this  integration  for  evaluating  a  and  b  was 

n  n 

carried  out  numerically. 

Previous  to  numerical  analysis,  if  the  expressions  for  a  and  b  are 

n  n 

evaluated  using  product  of  variable  integration  techniques,  chen  the  resulting 
expressions  become 

&n  =(2n+fj®TT®~  /  (iylj’ru<S»  m 


v, - H - 

n  "  (2n+l)sir3 


j  cp"(u)  sin 
0 


iTU  du 


where  cp"(u)  is  the  second  derivative  of  cp(u).  By  expressing  cp"(u)  by  a  finite 
difference  equation  where 


(70) 


cp"(u)  = 


(Au) 


a  (2Y.  -  Y.  -Y.  ) 

3  l  l-i  1+1 


and  numerically  evaluating  the  integrals  for  the  pulse  being  split  into  M- 
equally  spaced  intervals,  the  expressions  for  the  Fourier  series  coefficients 
are  given  by 


a 

n 


M 

2,jM  y 

(2n+l)s  A- 
'  '  i  =  0 


(2Y.  -Y. 
l  1-1 


-  Y.  )  cos 
1+1 


(2n+l)  i  t r 
j  M 


b 

n 


2,jM 

(2n+l)3  tt55 


-  Y.  -  Y.  ,  )  sin 
1-1  i+i 


/ 


( 2n+l)  i  tt 
jM 


(71) 


(72) 


Computer  programs  were  written  for  evaluating  Equations  63,  71, 
and  72  using  an  IBM  7044/7094  Direct  Couple  System  located  at  Wright- 
Patterson  Air  Force  Base  for  a  half  sine  type  pulse  as  indicated  in  Section  V. 
The  programming  was  accomplished  using  the  dimensionless  variables, 

r 

o  ,  and  B  in  determining  the  sinkage,  Z,  and  these  results  were  utilized 
V 

s  d 

in  the  development  of  the  flotation  parameter  nomograph  (see  Figure  17).  A 
complete  listing  for  the  computer  program  is  given  below. 


Table  XIII.  Dimensionless  Parameters  -  Uniform  Loading  Case 


a 

o 

“fi 

Fx 

F 

s 

cx 

k 

i 

0 

0.  21 

0.  0 

0.  99 

0.  000 

0.  62 

1.  00 

0.  25 

0.  205 

0.  03 

0.  96 

0.  141 

0.  62 

1.  00 

0.  50 

0.  20 

0.  065 

0.  94 

0.  305 

0.  63 

0.  97 

0.  75 

0.  185 

0.  09 

0.  87 

0.  42 

0.63 

0.  93 

1.0 

0.  16 

0.  115 

0.  75 

0.  54 

0.  64 

0.  88 

1.  5 

0.  115 

0.  150 

0.  54 

0.  71 

0.  61 

0.  64 

2.0 

0.  06 

0.  175 

0.  28 

0.  83 

0.  55 

0.  36 
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